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In this work we present laboratory measurements on the 
reduction of the threshold friction velocity necessary for 
lifting dust if the dust bed is illuminated. Insolation of 
a porous soil establishes a temperature gradient. At low 
ambient pressure this gradient leads to thermal creep gas 
flow within the soil. This flow leads to a sub-surface over¬ 
pressure which supports lift imposed by wind. The wind 
tunnel was run with Mojave Mars Simulant and air at 3, 6 
and 9mbar, to cover most of the pressure range at martian 
surface levels. Our hrst measurements imply that the in¬ 
solation of the martian surface can reduce the entrainment 
threshold velocity between 4% and 19% for the conditions 
sampled with our experiments. An insolation activated soil 
might therefore provide additional support for aeolian par¬ 
ticle transport at low wind speeds. 

1 Introduction 

Sand and dust are ubiquitous at the surface of Mars and 
in its atmosphere. Regional as well a s global dust storms 
are regularly observed fCantor l2007l l. In general, the at¬ 
mospheric dust is thought to play a major role for th e mar- 
tian climat e (Fed orova et al. l2014l : Montabone et al. l2014l : 
Kahre et al. l2013l l. On the groun d, eoli an features are wide 
spread. Recently, Bridges et al. |2012| reported that sand 
fluxes at the Nili Patera dune field on Mars are comparable 
to sand fluxes found on earth as the dunes move in almost 
steady state. This suggests that the underlying physics on 
Mars and Earth are similiar and that particle lift is related 
to saltation. However, at the low surface pressure moving 
particles by wind drag o nly is far from trivial. In earlier 
work Greeley et al. Il98d showed that most of the time the 
wind speed on Mars should be too low to move particles. 

Dust devils provide an additional view on parti cle lif t 
due to the pressure drop in their center (Balme et al. l200d 'l. 
The pressure difference between lower atmospheric pres¬ 
sure and higher sub-surface pressure with in the soil might 
support particle lift. Still, Neakrase et al. l20inl reported a 
significant discrepancy by a factor 3 between the observed 
mass fluxes and the estimation of dust entrainment at Gu¬ 
sev crater. Seen globally, the lifting related to dust storms 
migh t also dominate over lift by dust devils (Whelley et al. 
120081') . 

To compensate for this difference between observed dust 
flux and standard lifting capabilities of wind and dust dev¬ 
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ils, several further ideas were proposed. Gertain types of 
particle inventory like volcanic glass or aggregated parti¬ 
cles are more susceptible for pick up at lower wind speeds. 
The big, light-weight particles can roll over the surfaces 
which breaks the adhesi on before the p articl es are then 
lifted (de Vet et al. l2014l : Merrison et al. l2007l 'l. Gharging 
has also been discussed as a m echan ism to explain the high 
entrainment rates (Kok et al. l200fil l. If a number of parti¬ 
cles are ejected at low wind speed the threshold wind speed 
will also be lowered due to the reimpacting grains. Steady 
state saltation might then occur at lower wind speeds (Kok 
1201011 . For a revie w on sand and dust entrainment we refer 
to Merrison I 201 I 

It is somewhat counter intuitively, that owing to the low 
pressure "new" lifting mechanisms can be strong on Mars 
which are negligible at the 1 bar atmosphere on Earth. The 
mere insolation of a dust bed and related temperature gra¬ 
dients within the porous soil lead to str ong l ifting forces. 
This was already shown in Wurm et al. l2008l . though the 
underlying model changed somewhat (see below). In any 
case the insolation will "activate" the top soil layer. Ther¬ 
mal creep gas flow within the soil results from the temper¬ 
ature distribution and leads to sub-surface overpressure 
and therefore, as in the case suggested for dust devils, to 
an additional lift. This overpressure reduces the erosion 
threshold wind speed and enhances the entrainment rate. 
At specific geologic sites additional ga s flow s through the 
porou s ground occur (Antoine et al. I2OIII: Lopez et al. 
I 2 OI 2 I L This might also add a lifting force. 

Here, we report on wind tunnel experiments where we 
measure for the first time the reduction of the threshold 
friction velocity in dependence of an insolating light flux. 
We start by combining the illumination based lift with an 
aeolian lift model. We then describe the wind tunnel, the 
particle sample and experiments on particle lift. The data 
are used to give a measurement of the boundary layer pro¬ 
file and the likelihood for detachment in a light flux / wind 
velocity parameter set. We apply the lift model to this data 
set to obtain a threshold friction velocity depending light 
flux. This model fit is then scaled to martian conditions. It 
is used to quantify the reduction in friction velocity needed 
to lift particles on Mars due to insolation. 

1.1 Insolation Activated Layer 

In recent years a number of papers showed that at low 
pressure (hPa) particles can be ejected from a dust sam¬ 
ple b;j;_mere illuminati on wi th visible radi ation (Wurm et 
al. l200fil : Kocifaj et al. 1201 ll : Kelling et al. I2OIII : de Beule 
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et al. l2014h . This mechanism works independently of any 
wind if the illumination is strong enough, de Beule et al. 
I 2 OI 4 I found in microgravity experiments that this is caused 
by significant gas flow within the soil. This fact was not 
recognized before as thermal convection obscured this flow 
under gravity and only the microgravity experiments re¬ 
vealed the sub-surface gas flow. This flow within the dust 
bed is driven by thermal creep which pumps gas from a 
cold to a warm region if capillaries are present in a size 
comparable to the mean free path of the gas molecules. 
Mean free paths are similar to pore sizes in a dust bed at 
hPa pressure, hence we expect a large gas flow on Mars 
but an insignificant contribution under Earth’s high atmo¬ 
spheric pressure. 

A simple consideration of thermal cree p und er steady 
state conditions was proposed by Knudsen [l909|. By look¬ 
ing at the number of particles crossing the cross section 
of a capillary between two chambers with different tem¬ 
peratures Ti, T 2 a relationship between temperature and 
pressure can be deduced. The number of particles crossing 
the section h is given by 


pv 
noc — 


( 1 ) 


where v is the thermal velocity which is proportional to 
^/T. In equilibrium the gas flow from both sides is the 
same and 
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If the capillary is small, the overpressure exists as the gas 
only creeps along the wall and the pressure-driven back- 
flow cannot pass the capillary easily, but as the capillary 
becomes larger the flow due to the pressure difference can¬ 
cels this effect. 

Based on flow in small tubes the mass flo w rate M in a 
dust bed ca n be approximated as (Sone et al. Ilflflfll : Muntz 
et al. l2002l’l 
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with Pavg and Tavg as the average pressure and tempera¬ 
ture, P as porosity, A as surface area, fcs as Boltzmann 
constant, m as the molecular mass of the gas, Lr and 
as the radius and length of the capillaries describing the 
dust bed. AT and Ap are the temperature and pressure 
differences over the capillaries, respectively. The coeffi¬ 
cients Qp and Qt depend on the Knudsen number and 
describe the pressure driven (back) flow and the flow by 
thermal creep. The Knudsen number (Kn = A/d) is the 
ratio between the mean free path A of the gas molecules 
and a characteristic dimension d, which is Lr in this case . 
The coefhcients, Q, were adopted from Muntz et al. |2002| . 
We use an analytic relation fitted to these data (Koster, 
private communication) 


Q (Kn) = (Kn) = ,, , , , (4) 
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large scale net gas flow is present, which - again - trans¬ 
lates to the necessity of low pressure for the effect to be 
significant. A typical temperature profile of an illuminated 
dust bed is sketched in figure[T] Absolute values depend on 
the insolation light flux but calculated tempera ture p rofiles 
always show a similar structure (Kocifaj et al. 1201 il l . 

On one side the gradient below the surface implies that 
thermal creep pumping drives a gas flow upwards. On the 
other side the uniform temperature over the top layer im¬ 
plies that gas flow can only be maintained if a pressure 
difference is present as given by Darcy’s law (also included 
in eq. [3] as pressure term). de Beule et al. 1201, ’ll call this 
pressure supported top layer an insolation activated layer. 
While we do not use their results they found the activated 
layer thickness to be between 100 and 200 /im in their ex¬ 
periments though this will vary with the environment or 
experimental settings and parameters. In any case their 
experiments confirm the existence of a sub-surface over¬ 
pressure. Balancing the pressure (AT = 0) and thermal 
creep (AP = 0) driven mass flows leads to an overpressure 
of 
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with the thickness of the layer with the temperature gra¬ 
dient Lxi and of the uniform temperature layer La: 2 . 

We parameterize the associated lifting force Fai on 
grains within the activated layer as 


Fal = Call 
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where Cai is a constant which characterizes the thermal 
and optical properties of the dust bed. The other param¬ 
eters are kept explicitly for later scaling to martian condi¬ 
tions. I is the light flux generating AT. As a first order 
approximation we consider both to be proportional to each 
other. 


1.2 Aeolian Lift Model 

Aeolian entrainment is complex and cannot be calculated 
from first principles. Therefore the existing models of en¬ 
trainment consist of semi-empirical laws formulated for 
Earth. To apply these on Mars the y are rescaled to mar- 
tian conditions (eg. Greeley et al. llQSd : Merrison et al. 
l2007l f . The force which a free particle with radius r expe¬ 
riences in a gas flow due to dynamic pressure is 

Fd = \Cdpu\\ (7) 

where u is the flow velocity, p is the gas density, and Co 
is the drag coefhcient, which depends on the flow charac¬ 
teristics. The force on a free particle would act along the 
flow direction. 

At the soil / atmosphere interface the flow and its in¬ 
teraction with the surface leads to a gas drag lifting force 
opposing gravity. This force is supposed to have the same 
dependencies as the drag on a free particle, namely 

Fl = ClPuW (8) 


with ai = 0.5007 ± 0.0001, 02 = 0.990 ± 0.006, 03 = 
0.743 ± 0.002, 04 = 0.222 ± 0.004, 05 = 1.241 ± 0.004. 
For "large" pores the pressure back flow dominates and no 
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It is characterized by a lifting coefficient Cl which holds all 
unknowns and has to be a fraction of the drag coefficient 
Cd. 































Figure 1: Schematics of the temperature in an illuminated dust bed. The top layer of the dust is heated by the 
illumination. Gas from the bottom (or non illuminated regions) moves due to thermal creep. This leads to 
an overpressure just under the surface as in the topmost layer the gas can only escape by Darcy flow. 


If the flow is not uniform across the particle diameter, 
the flow velocity gradient will also generate a torque, which 
can help to reduce the adhesion and can initiate the rolling 
motion of larger grains. The force related to this torque 
can be described by 

Ft = Ct pu^ , ( 9 ) 

noting that Ct is not dimensionless but has the unit m“^ 
now. Cd is well known for spheres, the lift and torque 
coefficients, however, are not well constrained. 

1.2.1 Surface Flow 

The gas flow at the surface results in a the shear stress r. 
Related to this a shear velocity is defined as 

( 10 ) 

The flow in the vicinity of the surface follows a logarithmic 
profile 

u{z) = — In— (11) 

K, Zq 

with the Karman constant k = 0.41 and the surface rough¬ 
ness Zq. Gas drag depends on the friction velocity. 

1.2.2 Slip Flow Regime 

At large Knudsen numbers molecules passing a particle 
closely keep a ballistic trajectory (slip flow). At lower 
Knudsen numbers they are slowed by other molecules de¬ 
flected from the particle surface but the velocity compo¬ 
nent parallel to the surface does not vanish completely. For 
this regime the drag can be quantified by the Stokes law 
with an additional slip flow correction factor depending on 
the Knudsen number 

Fo = —67r/iru/C(Kn). (12) 

with fi as viscosity. The Cunningham correction factor C 
is a semi-empirical function of the Knudsen number 

C(Kn) = 1-f + 52 exp Kn (13) 

with gi = 1 . 207,52 = 0.44 and 53 = 0.78 (Rader HoflOl l. 



1.3 Combining Aeolian Lift with the Activated 
Layer 

We are interested in the dependency of the light flux 
needed to initiate dust motion at a given shear velocity. 
The erosion threshold at a given light flux is defined by 
the balance between lifting forces Fl and Fai on the one 
side and the opposing forces gravity Fq and cohesion Fc 
or 

CL'^pr^ul + Cail^^^Fa + Fc. (14) 

2 1 Qp 

For a given grain diameter (not varied here) the lifting 
related to torque can effectively be included in Cl- We 
therefore did not explicitly trace this term further on but 
this would need to be included if particle variations would 
be considered. As critical light flux Frit depending on the 
velocity m* we therefore get 


Fg+Fc -CL^pr^ul 

Cal^Q 


(15) 


It has to be noted, that Q depends on the Knudsen num¬ 
ber, and Cl is a function of the flow parameters. Two 
choices are possible: the usual one is taking Cl to be a 
power law function of the Reynolds number. On the other 
hand one can use the slip flow corrected equation. The 
drag coefficient Go is already inserted into the slip flow 
equation, only the proportionality constant between drag 
and lift is still unknown. 


2 Wind Tunnel Setup 
2.1 Wind Tunnel Capabilities 

To study the influence of insolation activated lift on the 
threshold wind speed we constructed a new w ind t unnel 
based on an earlier set up by Paraskov et al. I 2 OO 6 I . The 
wind tunnel has a total length of 18 m divided in sev¬ 
eral horizontal and vertical parts, including one horizontal 
and two vertical experimenting sections (upward, down¬ 
ward and horizontal wind), sketched in figure[21 It consists 
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of stainless steel tubing with an inner diameter of 32 cm. 
The channel can be evacuated to and work at pressures 
between 10 Pa and 1000 Pa. The wind is provided by a 
roots pump circulating up to 8m®/s. This translates to a 
wind speed averaged over the cross section up to 100 m/s. 

2.2 Experimental Setup 

The experiments were carried out in the horizontal 2.4 m 
long experimenting section. 

A sample holder with a simple wing shape is suspended 
in the middle of the channel by several threads. A dust 
sample with a depth of 2.7 cm and a square surface of 
92 cm^ can be prepared and inserted. The dust sample 
is placed in the center of the channel as sketched in figure 

m 

A cw infrared laser (1064 nm) at varying light flux illumi¬ 
nated the dust bed. To visualize the (airborne) particles an 
additional red laser curtain was implemented. Dust parti¬ 
cles are detected by a long term recording camera running 
at a rate of 100 frames per second, with maximal exposure 
time (9994 ps). The heating of the dust bed is dominated 
by the IR laser while the heating due to the laser curtain is 
insignificant. As its scattered light would overexpose the 
images a line filter in front of the camera removes the IR 
and only the red laser curtain is visible. The laser cur¬ 
tain was pulsed with 3 short pulses of 164 ps duration with 
176 ps delay in between followed by a 2.45 ms break, so ev¬ 
ery particle shows up in the camera pictures as a distinct 
pattern from which the speed and acceleration could be 
calculated. 

In this work we only distinguish whether lift occurs in a 
given experiment. A measurement of lifting rates is out¬ 
side the scopes of this paper. The infrared laser allows for 
an intensity much higher than solar, which is not meant 
to simulate insolation on Mars but improves the determi¬ 
nation of the dependence of the threshold light flux on 
friction velocity at martian values. 

2.3 Sample Preparation 

As dust sam ple m aterial we used Mojave Mars Simulant 
(Peters et al. l2008l j. It consists of crushed basaltic rock. It 
is more angular than the aeolian weathered dust expected 
on Mars. As it is not weathered by the influence of aqueous 
alteration it has a low volatile content of only up to 4 %. It 
was additionally treated by heating it to 400 °C for 6 hours, 
to reduce the volatile content further. The size distribution 
of particles in the dust bed was measured with a commer¬ 
cial particle sizer (Mal vern M astersizer 3000) and is shown 
in figure[3l Peters et al. l2008l consider this dust to be repre¬ 
sentative for the mechanical properties of the fine grained 
material on Mars, as the inventory of sub micron particles 
is still included in the sample. It has also spectral charac¬ 
teristics and material comp osition which are close to those 
of Mars (Peters et al. l2008l ~l. Therefore, it was deemed as 
a suitable sample material for a first study here. 

The dust was sieved into the sample holder to achieve a 
loose settling of the grains and prevent compaction which 
may alter the physical properties of the surface. In this 
process some very loosely bound aggregates or grains may 
reside at the surface, and those particle are lifted easily. 



Figure 3: Volume size distribution of the dust samples. 

Peak of the distribution is at 27 pm and 50 % 
of the dust is smaller than 56 pm. 

We did not intend to study the threshold speeds for ero¬ 
sion of some particles which are exceptionally susceptible 
to pick up but were interested in the more constant en¬ 
trainment thresholds. As procedure to remove these parti¬ 
cles and reproduce similar conditions for all experimental 
runs the sample was passivated before the measurement. 
This was achieved by a >30 s gas flow with a higher speed 
than used in the measurement. This removed all particles 
that are only loosely bound to the surface and allows the 
determination of a threshold for continuous entrainment 
afterwards. This is not necessarily the most natural soil 
expected on Mars where atmospheric dust on the order of 
1 fim might settle on top of a coarser bed. However, it is 
a well defined way of preparing the dust bed and warrant 
similar conditions for each experiment run. The ambient 
pressure was set to 320 Pa, 630 Pa and 930 Pa, with uncer¬ 
tainties of 5 % during the experiments, and air was used as 
working gas. 

In each experimental run a given sample is observed at 
different flow velocities and light fluxes. First the lowest 
light flux is set and the wind speed is increased. Then 
speed is decreased again, the light flux is increased and so 
on. This gives a matrix which samples the 2d space of wind 
speed and light flux. Each measurement consists of 2750 
frames (27.5 seconds of data acquisition at 100 fps). These 
frames are processed and overlayed to produce one image 
for each measured sequence which shows all entrained par¬ 
ticle - if present. This way a decision can be taken for each 
image sequence if particles are picked up within the laser 
curtain or not. An example is given in figured] 

3 Results 

The data are shown in table [T] Here, the flow is character¬ 
ized in terms of friction velocity as described in the next 
section. 

3.1 Boundary Layer Profile 

To determine the boundary layer profile, the velocity of 
the entrained grains was measured in dependence of the 
height. We used one data set at the highest wind velocity 
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Figure 2: Schematics of the wind tunnel. The wind tunnel can be evacuated by a vacuum pump. The roots pump 
generates the air flow. As an inset the experiment section is shown. In this section the sample is placed in a 
suspended sample holder and is irradiated by an infrared laser. The IR-radiation is blocked from the camera 
(recording with 100 fps) by a red line-filter. Particles are visualized by additional illumination with a pulsed 
red laser curtain. 



Figure 4: Example of a processed sequence of images showing lifted particles. 
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Table 1: Experiment results for three used pressures. The number of experiments (out of 6), where particle lift occurred, 
is shown depending on light flux (first column) and friction velocity. The thresholds were determined where at 
least two measurements showed erosion (or the previous measurement already showed erosion). The thresholds 
are shown in bold. Towards lower pressure the rates become smaller and the threshold is less visible, therefore 
the data at 205 Pa was not used for fitting but is shown for completeness. 
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Figure 5: Velocity over height for grains detached at 
largest wind speed used at 320 Pa. A fit to the 
data according to equation IIII is shown. 

used at 320 Pa. Here a large number of particles could be 
detected to allow a flow characterization. The last visi¬ 
ble occurrence of each particle was used as a data point. 
There the particles were already well coupled to the gas 
flow and hence trace its velocity. The boundary layer flow 
is shown in figure O The data were fitted to equation 1111 
to obtain surface roughness and friction velocity. The fric¬ 
tion velocities at other settings was calculated assuming 
that M, oc Umean- The average flow velocity in the channel 
Umean is known for the pump system. 


Figure 6 : Model fits to the entrainment thresholds ob¬ 
tained at 930 Pa, 630 Pa and 320 Pa. The data at 
205 Pa was not used for fitting. The Newtonian 
drag fit is shown as solid lines and the slip flow 
fit as dashed lines. Both models have a similar 
fit quality. For details we refer to the text. 

four parameters obtained are: <5; at = 0.59 ± 0.24 kWs^/m^, 
7 = —0.83 ± 0.07, dN = 32 ± 11 pm, and Fgc,n = 1.20 ± 
0.46 kW/m^. The parameter dN =32pm is consistent with 
the particle size. 7 gives the dependence on the Reynolds 
number which is oc with -l<x<0. Therefore, the value 
for 7 is plausible. 

For the slip flow the fit function is given by 


3.2 Detachment 


/crit ^ + Fbc.c^ , (17) 


All measurements are summarized in table [T] This shows 
the number of runs where lifted particles were detected de¬ 
pending on wind speed and light flux, 6 runs were made 
for every combination thereof. It is clearly visible that the 
threshold speed is reduced with increasing illumination and 
decreases with decreasing pressure. At 205 Pa the deter¬ 
mination of the threshold is ambiguous and this data was 
not used for any fitting but is shown for completeness. 

The data have to match the dependency of equation 1 151 
As QtIQp depends on the Knudsen number, the mean 
free path A was calculated but the particle size d was left 
as a fit parameter. For the lifting coefficient C7n(Re) two 
approaches were used and compared: The ansatz in the 
literature is using mostly Newtonian drag and assuming 
a power law function for Cl- A second ansatz was also 
tested. For the given particle size and velocities the ap¬ 
propriate drag should be described by the slip flow regime, 
therefore this dependency was also used. 

For the Newtonian drag a power law Cn(Re) = const • 
Re^ was assumed, where the constant is combined with 
the other parameters in front of uj as new parameter ajv- 
Two constants m*,o = Im/s and po=600Pa were inserted 
to avoid confusion with the units. The parameters {Fq -I- 
Fc)/CaiT were replaced by one fitting parameter Fgc,n- 
In total this gives the function 

hru {n„p) = fQ{^) “• + > ( 16 ) 

which was fitted to the data obtained for all pressures. The 
fit to our experimental data can be seen in figure [5] and the 


with the three fit parameters being ac = 0.0086 ± 
0.0026 kWs/m'‘,dc = 56 ± 7pm, and Fqc.c = 0.70 ± 
O.lkW/m^ 


3.3 Scaling to Mars 


The laboratory results have to be scaled to martian condi¬ 
tions to account for differences in insolation, gravity, tem¬ 
perature and gas composition. The differences between 
laboratory conditions and Mars are summarized in table 

m 

Every parameter of equations 1161 and HZ] can be recaled 
by multiplying it by a scaling factor /. For the New¬ 
tonian drag we use the following: QtIQp is a function 
of the Knudsen number. The Knudsen number scales 
with the mean free path A = (an)~^ where sigma is the 
molecular cross section and n = p/ki,T the number den¬ 
sity. Therefore is can be calculated that A oc T/cr and 
f\ = 0.923. The insolation induced lifting force scales 
with 1/T, and fuft = 1.35. The drag force scales with the 
density p oc p- mjT where m is the molecule mass and the 
scaling factor takes the value fp = 2.06. The drag coeffi¬ 
cient scales with the Reynolds number, which again scales 
with the density. These factors can be inserted to rescale 
the fit, giving 


(18) 

Only /fgc difficult to obtain as the degeneracy between 
gravity and cohesion can not be disentangled from this 
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Table 2: Comparison of relevant laboratory and martian parameters, the ratio is the value on Earth divided by the 
martian value. 


Quantity 

Unit 

Symbol 

Earth 

Mars 

ratio [%] 

gravity 

[m/s^] 

g 

9.81 

3.69 

37.6 

pressure 

[hPa] 

P 

1013 

6 

0.6 

temperature 

[K] 

T 

295 

218 

73.9 

major Atmospheric component 



Na 

CO 2 

- 

molecular mass 

[u] 

m 

28.96 

44.01 

152 

molecular cross section 

m^] 

<7 

11.34 

9,08 

80.1 

max insolation 

[kW/m^] 

imax 

1.412 

0.717 

50 


dataset. Therefore the two limiting cases will be treated: 
If Fgc is dominated by gravity /f^c = 1/3 as gravity on 
Mars is only 1/3 compared to Earth. If cohesion dominates 
fpcc = 1 as cohesion stays the same. The real value will 
be somewhere in between. The Cunningham fit function 
and the other parameters can be scaled accordingly if slip 
flow is considered instead of Newtonian flow. 

For the maximum insolation on Mars of 717 W/m^ we 
get threshold friction velocities compiled in table O Com¬ 
pared to the calculated threshold friction velocity with no 
insolation the reduction lies between 4% and 19%, de¬ 
pending on wheter the holding force is cohesion or gravity. 
The rescaled fit is shown in figure [T] 

Under our laboratory conditions Newtonian flow and 
corrected slip flow models differ only significantly at very 
low pressure. Due to the limited capability of particle lift¬ 
ing though, the measurements do not allow to distinguish 
between both. Scaled to Mars the differences get some¬ 
what larger. However, it should be noted that the relative 
reduction is rather robust independent of the drag law con¬ 
sidered. 

4 Discussion 

The manually prepared dust beds and artificial illumina¬ 
tion can provide a model of a planetary surface, but not 
all parameters can be mimicked. The Reynolds number 
at 600 Pa and 30 m/s mean velocity for our wind chan¬ 
nel is 3200. The 30 m/s correspond to a friction velocity 
of 1.86 m/s. This means that the flow is slightly above 
the critical Reynolds number for some measurements - 
namly the highest velocity at 320 Pa and the two highest 
at 205 Pa - and hence should be prone to become turbu¬ 
lent. At lower speeds for higher light flux it drops below 
the critical Reynolds number though. We do not expect 
that the transition influences the systematics of the mea¬ 
surements, as the measurement at 205 Pa was not taken 
into account. 

The IR light source illuminated an area of 1.77 cm^ and 
was on during the experiment run lasting about half an 
hour. This implies a spatially confined and time dependent 
temperature variation over the course of each experiment. 
Gas flow and pressure build up below the surface might 
therefore vary over time. To what extend this will affect 
the data is unclear. However, as the expected dependence 
fits the data rather well we do not expect the systematics 
to change strongly. 

Care has been taken to increase the light flux during a 
measurement, as decreasing the light flux can otherwise re¬ 


sult in particle ejection. This is a re lated lifting mechanism 
e.g. studied in Kelling et al. l201ll . It is based on a differ¬ 
ent temperature profile within the soil developing after a 
decrease of the illuminating flux. This might be of high 
importance for particle lift e.g. in dust devils providing 
shadowing on short time scales or a shadow cast by a sur¬ 
face feature moving over the ground. However, this is not 
the topic of this work, but shall be examined in the future. 
Here we are interested in the supporting lift for a continu¬ 
ous illumination in face of a wind flow. By increasing the 
light flux during the experiments we avoid the occurrence 
of such particle lift. If a dust storm sets in, its opacity will 
hinder the insolation, for a sharp cut in optical density this 
might add to the lifting. But as saltation was not present 
at this experiment - as the saltation length is larger than 
the experiment section -, the evolution of a dust storm can 
not be predicted here. 

There is some variability between the samples, where 
some samples are more readily eroded than others but this 
has to be expected for dust samples prepared manually. 
This raises the question how a natural dust surface would 
behave if small dust particle in the micrometer range would 
settle on top of a dust bed. As this influences the gas flow 
and pressure profile in the dust bed it might have some 
influence on particle lifting. The permeability of the top 
layer gets smaller which requires a higher pressure to force 
the gas flow through. This would ease particle lift. This 
has to be studied in future experiments. 

Note that we observe no signs of electrical charging ef¬ 
fects here and do not consider it to play a role in the con¬ 
text of this study, as the amount of charge will not be 
influenced by the infrared radiation, and cannot provide a 
relative reduction in respect to the illumination. 

Direct predictions for the dust entrainment activity is 
not straight forward, as a multitude of properties have to 
be taken into account. The highest intensity available is at 
perihelion at sub-solar latitude. But temperature and gas 
density play major roles. Therefore, dust activities should 
be tied to certain regions on Mars, which is consistent with 
obser vation s. This correlation was already found by Wurm 
et al. l2008l . who investigated the illumination effect alone. 

As the top layer insolation activation depends on nu¬ 
merous parameters the absolute wind speed reduction is 
variable. Pore sizes and particle size distribution are only 
two examples to be varied. We do not deem it possible to 
quantify the effects of all variations here but several options 
would increase the activation effect. Further experiments 
would have to be carried out to quantify the influence of 
the different parameters. 














Figure 7: Threshold conditions rescaled to martian environment; solid lines: cohesion dominates over gravity, dashed 
lines: gravity dominates over cohesion; note the different scales on the x-axis. 
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Table 3: Threshold friction velocities under different conditions on Mars 

Newton 



lab values 

fPcc = 1 

fFac = 1/3 

p 


jlmax 

red. 



red. 



red. 

[Pa] 

[m/s] 

[m/s] 

[%] 

[m/s] 

[m/s] 

[%] 

[m/s] 

[m/s] 

[%] 

205 

2.14 

2.07 

3.6 

1.94 

1.85 

4.6 

0.75 

0.65 

14.1 

320 

2.02 

1.93 

4.2 

1.81 

1.73 

5.2 

0.78 

0.59 

16.1 

630 

1.83 

1.75 

4.4 

1.66 

1.57 

5.3 

0.64 

0.54 

16.3 

930 

1.74 

1.67 

4.0 

1.57 

1.50 

4.8 

0.61 

0.52 

14.9 

Slip Flow 


lab values 

II 

0 

fPac = 1/3 

P 



red. 

o 

11 

* 


red. 



red. 

[Pa] 

[m/s] 

[m/s] 

[%] 

[m/s] 

[m/sj 

[%] 

[m/s] 

[m/s] 

[%] 

205 

2.60 

2.47 

4.9 

2.50 

2.35 

6.0 

0.83 

0.68 

18.5 

320 

2.16 

2.05 

5.0 

2.10 

1.97 

6.1 

0.70 

0.57 

18.8 

630 

1.81 

1.73 

4.4 

1.78 

1.67 

5.2 

0.59 

0.50 

15.9 

930 

1.70 

1.63 

3.7 

0.56 

0.49 

4.3 

0.56 

0.49 

13.2 


However, these first experiments imply that insolation 
activation can put significant tension on a top soil layer 
under the low pressure conditions of Mars and enhance par¬ 
ticle lifting by wind or within dust devils. In this first sys¬ 
tematic analysis of insolation aided wind erosion we found 
a reduction of the threshold speed of about 4% to 19 % 
for martian conditions. As it is unlikely that we acciden¬ 
tally probed the most susceptible conditions, this number 
is expected to be higher for some parameter combinations. 
Therefore, the insolation activated layer might strongly en¬ 
hance particle lifting at lower wind speeds. 

5 Conclusion 

In this work we present laboratory measurements on the 
reduction of the threshold friction velocity necessary for 
lifting dust if the dust bed is illuminated. Insolation driven 
thermal creep can activate a dust layer by adding a sub¬ 
surface overpressure. This supports lift imposed by wind. 
Our first measurements imply that the insolation of the 
martian surface (with 717 W/m^) can reduce the entrain¬ 
ment threshold velocity between 4 % and 19 % for the con¬ 
ditions sampled by the experiments. An insolation acti¬ 
vated soil might therefore provide additional support for 
aeolian particle transport at low wind speeds. 

The reduction of up to 19 % is one step to bridge the 
gap between observed thresholds and theoretical predic¬ 
tions though it might not explain all discrepancies. It still 
has to be investigated how this mechanism works out with 
sand sized particles or a broad size distribution of parti¬ 
cles. Chaning pore size and adding particles susceptible to 
lift might be promising in reducing the gap further. 
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